A DUALITY METHOD IN PREDICTION THEORY OF 
MULTIVARIATE STATIONARY SEQUENCES 



BY MICHAEL FRANK AND LUTZ KLOTZ OF LEIPZIG 



Kurztitel: Prediction theory of stationary sequences. 

2000 Mathematics Subject Classification. Primary 60G25, 60G10; Secondary 42A10 



Keywords and phrases. Multivariate stationary sequence, prediction problem, non-nega- 
tive Hermitian matrix-valued weight function, trigonometric approximation. 



Abstract. Let W be an integrable positive Hermitian q x q- matrix valued function on 
the dual group of a discrete abelian group G such that W~ x is integrable. Generalizing 
results of T. NAKAZI and of A. G. MlAMEE [|N| and M. POURAHMADI ||MiP|| for q = 1 we 



establish a correspondencebetween trigonometric approximation problems in L 2 (W) and 
certain approximation problems in L 2 {W~ l ). The result is applied to prediction problems 
for g-variate stationary processes over G, in particular, to the case G = Z. 
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1. Introduction 

In 1984 T. Nakazi [N) introduced a new idea into prediction theory of univariate 
weakly stationary sequences. Under the additional assumption that such a sequence has an 
absolutely continuous spectral measure and its spectral density w is such that w' 1 exists 
and is integrable he related approximation problems in L 2 (w) to certain approximation 
problems in L 2 (ty _1 ). His method opened a way for him to give an elegant proof of 
Szego's infimum formula and to obtain a partial solution of a certain prediction problem 
which will be called Nakazi's prediction problem in the present paper. A. G. Miamee 
and M. Pourahmadi [ MiPjl pointed out that the essence of Nakazi's method consists in 



a certain duality between the Hilbert spaces L (w) and L (w ). This way they found 
a unified approach to several 'classical' prediction problems, and they generalized some 
of Nakazi's results, partially even to more general harmonizable stable sequences. The 



papers ||Mi] , |CMiP|| contain further completions of these results. 

The aim of the present paper is the application of Nakazi's duality method to g-variate 
(q G N, the set of positive integers) weakly stationary processes over discrete abelian 
groups G. Under the assumption that the inverse W~ x of the q x g-matrix valued spectral 
density W of such a process is integrable we establish duality relations between the left 
Hilbert modules L 2 (W) and L 2 (PU~ 1 ). Section three of our paper contains the general 
results. The further sections are devoted to applications. 

In Section four we obtain some well-known prediction results by quite different proof 
methods. We compute the one-point interpolation error matrix and derive Yaglom's inter- 
polation recipe under the mentioned assumption on W. Moreover, Section four contains 
our generalization of Nakazi's proof of Szego's infimum formula to the multivariate situ- 
ation. This way we obtain a special case of a result due to V. N. Zasukhin [0] as well 
as to H. HELSON and D. LOWDENSLAGER llHelLL Thm. 81. 



If G = Z is the abelian group of integers and the index set S of the known values 
consists of all negative integers and the set of integers {1,2, ...,n} for some n G N then 
the arising prediction problem is called Nakazi's prediction problem. In Section five we 
solve this problem under the additional assumption that log det W is integrable. We 
will see that in case W~ x is integrable the result is an easy consequence of our duality 
results. The more general case in which merely log det W is integrable can be solved by 
approximation procedures. Section five also contains some straightforward multivariate 
generalizations of univariate results of [MiPj . 



2. Preliminaries 

Let G be a discrete abelian group with neutral element 0, G* be its dual group, and 
A be the normalized Haar measure on G*, i.e. A(G*) = 1. Relations between measurable 
functions on G* are to be understood as relations which hold true almost everywhere 
(abbreviated to "a.e.") with respect to (abbreviated to "w.r.t.") A. Integration is always 
ment to be done over G*. For a subset S of G set S° = S U {0} and S c = G \ S°. For 
q G N denote by M g the algebra of all complex- valued q x g-matrices. The zero matrix 
and the identity matrix of M g are denoted by and /, respectively. For a non-empty 
subset S of G let T(S) be the left M g -module of all M g -valued trigonometric polynomials 
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with frequencies of S. Functions of T({0}), i.e. constants, and elements of M g will be 
identified and denoted by the same symbols. 

If A G M y then the symbol A* stands for the adjoint of A, det A for its determinant, tr A 
for its normalized trace and \A\ for its normalized euclidean norm, i.e. \A\ 2 = tr(AA*). 
If A is regular then A -1 denotes its inverse. The set of all Hermitian matrices will 
be equipped with Loewner's semi-ordering. In particular, a maximum or minimum of 
a subset of the set of Hermitian q x g-matrices is to be understood w.r.t. that semi- 
ordering. The cone of positive Hermitian q x q- matrices is denoted by M>. For A G M> 
we denote its unique positive square root by A 1 / 2 . Finally, let L 1 and L 2 be the linear 
spaces of (equivalence classes of) M 9 - valued functions on G* that are integrable or square- 
integrable, respectively. 

Let W 9 (G*) be the set of M> -valued functions of L 1 and W 9 (G*) = {W G W g (G*) : 
W~ l G L 1 }. For W G >%(G*) the symbol L 2 (W) denotes the left Hilbert M 9 -module of 
(equivalence classes of) M 9 -valued functions F on G* such that tr(FWF*) is integrable 
w.r.t. A. The M^-valued inner product on L 2 (W) is defined by 

(F, G)= J FWG* dX , 

the corresponding scalar product by 

(F, G) = tr(F, G),F,Ge L 2 (W). 

The derived Hilbert space norm on L 2 (W) is denoted by || ■ ||. For a more general con- 
struction see [[Ros|l . 



We recall that if M. is a (closed) submodule of L (W) and F G L (W) there exists 
a unique function Fm £ -M such that (F — Fm,F — Fm) is the minimum of the set 
{(F — G, F — G) : Gg M.}. The element Fm is the result of the orthogonal projection of 
F onto M. w.r.t. both (-, ■) and (■, ■), i.e. (F — Fm, G) — for any G G M.. The orthogonal 
complement of M is equal to the submodule {H G L 2 (W) : (H, G) = for all G G M}. 
The function (F — Fj^) is the result of the orthogonal projection of F onto the orthogonal 
complement of M.. A detailed study of the geometry of Hilbert M 9 - modules can be found 
in ||GH|| , for applications to prediction theory cf. [ WMa |. 



Considering L 2 (W) as the spectral domain of a g-variate weakly stationary process 
over G we can formulate linear prediction problems of such a process as trigonometric 
aopproximation problems in L 2 {W). If S C G let M.{S) be the closure of T(S) in L 2 (W) 
and Py be the orthogonal projection onto M.{S). If S C G \ {0} the determination of 
Pgl as well as of the prediction error matrix 

(2.1) A s = (/ - P S I, I - P S I) = (/ - P S I, I) 

is of fundamental importance in the prediction theory of g-variate weakly stationary pro- 
cesses over G. This kind of determination problem is sometimes called the general pre- 
diction problem. There exist more or less complete solutions of it for some special choices 
for the set S. For an introduction to prediction theory of g-variate weakly stationary 
processes we refer to ||Koz} . 
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3. Duality relations 

Let W G W 3 (G*). Then along with L 2 (W) one can consider the left Hilbert Mq-module 
L 2 (W- 1 ). The symbols (•,•), (v), || • ||, M(S), P s and A 5 referring to L 2 (W) will be 
replaced by (•, -)^, (•, -)^, || • ||^, M.(S), Ps and Ag, respectively, for the denotation of the 
corresponding objects related to L 2 (W /_1 ). It turns out that a certain duality between 
L 2 {W) and L 2 (W / ~ 1 ) is helpful if we study the general prediction problem in L 2 {W). 

Lemma 3.1. Let W G W q (G*). Then the mapping 

F — >FW , F G L 2 (W) , 
is an isometric isomorphism of L 2 (W) onto L 2 (W~ r ). 

Proof: UF,Ge L 2 (W), we obviously have (F, G) = (FW, GW)». If H G L 2 (W- 1 ) 
then HW~ l G L 2 (iy) and H W~ l W = H. □ 

Another immediate consequence of the integrability of W~ l is the point of view on 
L 2 {W) and on L 2 (W~ 1 ) as on subspaces of L 1 . This can be seen from the inequality 



(3.1) J \F\ dX = J \FW 1/2 W~ 1/2 \ dX < J \FW l/2 \-\W~ 1/2 \dX 

r \ 1/2 / f \ 1/2 

J {FW^dXj (J {W-^dX 



< 

< oo , F G L 2 (W) 
If Ai is a submodule of L 2 (W) then the submodule 

M d - 



G G L 2 (W~ l ) : J FG* dX = , F G M j C L 2 ^" 1 ) 



is said to be its dual. Note, that Ai d is the orthogonal complement in L 2 (W 1 ) of the 
submodule MW = {FW : F G M}. For a subset S C G \ {0} the symbol J s denotes 
the orthogonal projection of / onto Al(5' ) c( C L 2 (iy _1 ). 

Theorem 3.2. Lei G VV,(G*). Then for subsets S C G \ {0} 

(3.2) P s 7 = / - (7 - J s , /)-!(/ - J^iy- 1 
and 

(3.3) A s = (/ - J s , J): 1 = (I — Is, I — IsV ■ 
In order to prove Theorem [372] we need the following lemma. 

Lemma 3.3. Let F G M(S°). Then F G M(S) if and only if f F dX = 0. 

Proof: If F G T(S') then f F dX = 0. If F G Al(S') we approximate it by elements 
of T(S), and by ( p.l|) we obtain f F dX = 0. If F G A^S* ) it can be written in the form 
F = A + G, where A G M and G e M(S). By the result just proved we have / F dX = A. 
Consequently, J F dX = implies F = G eM(S). □ 

Proof of Theorem Since J 5 G A< (S°) d we have JI s dX = f I S I* dX = 

and hence, (J — Is, W)^ — J(I — Is) dX = I. By a generalization of the Cauchy-Schwarz 
inequality the inequality {W, W)Z 1 = {I-Is, W)^{W, W)- l {I-I s , W)~ < (I-is, I-Is)~ 
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holds, cf. |3|, Thm. 1] . Therefore the right hand side of ( |3.2| ) is well-defined. Let us denote 
it by F s . Since (I - I s ) G M(S°)W we conclude F s G M(S°) by Lemma |3J]. Because 
JF s dX = I-(I- I s , I)Z\I ~ Is, J)~ = Lemma U implies 

(3.4) i^eA^S). 

Let F G .M(S). Then (I,FW)„ = J F* dX = by Lemma [T3J Therefore, the equality 
(/ - F S ,F) = (I- I S ,I)Z\I ~ Is,FW)„ = (I- Is,I)z\I,FW)^ = yields. By a 
combination with ( |3.4j ) we get Fs = PsI- To obtain formula ( |3.3j ) simply replace one 
term according to ( ^.2| ) in (|2.1|). □ 

From Theorem |3.2| we can derive a general result which seems not have immediate 
applications to prediction theory of multivariate stationary processes, however which can 
be considered as a certain generalization of the univariate case. 

Let Mg be the set of all gx g-matrices whose determinant equals to one. For S C G\{0} 
define 

(3.5) 5 S = inf{p - F\\ 2 : A G M.' q , F G . 

Lemma 3.4. Let W G W g (G*) and 5 6e a stz&sei o/G\ {0}. Then 

(3.6) 5 S = [det^)] 1 / 9 . 

Proof: We have the equality [det(J — F,I — F)] 1 ^ = inf{tr((J - F,I - F)B) : B G 
M> n M^)}, F G M{S), cf. flHol , Section 7.8, Problem 19]. Since B = A* A for some 
A G .M' we obtain [det(J — F, I — F)] 1 ^ = M{\\A - AF\\ 2 : A G Af}, F G M(S), and 
therefore 

[det(A s )] 1/9 = {det[min{(/ -F,I-F):Fe M(S)}]} 1/q 
= min{[det(J — F,I — F)} 1 ^ : F G M(S)} 
= inf{p - AFf : A G M; , F G 

= 

□ 

Remark 3.5. Since a matrix A G can be written in the form A = UB for a unitary 
matrix U and for a matrix F G (M> D M' g ) the value of 5s does not change if ^4 runs 
through (M> H M') on the right-hand side of ( |3.5|) . Moreover, if A exhausts the larger 
set of all q x g-matrices with | det A | = 1 we get the same result. 



Combining Theorem |3.2| and Lemma |3.4| we get the following assertion. 
Theorem 3.6. Let W G VV g (G*). Then for subsets S C G \ {0} we have the equality 
(3.7) 5 S = (inf{||A - G\\i : A G M' , G G M(S°) d })- 1 . 



Proof: In a similar way to that one taken in the proof of Lemma |3.4| we derive the 
equality 

M{\\A-G\\i : A G M; , GeM(S°) d } = {det[min{(I - G, I - G)^ : G e M(S°) d }]} 1/(1 

= {det(I-I s ,I-Is)~} 1/q - 
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By (|3.3|) and ( |3.6| ) the right-hand end of this equality equals to 5$ . □ 
Another duality relation can be derived by a simple adaptation of a chain of equalities 
mentioned in the proof of 0, Thm. 4.1] to the multivariate situation. 

Theorem 3.7. Let W G W q (G*). Then for subsets SCG\{0} one has 



(3.8) inf{||A - P S I\\ : A G M 9 , tr A = 1} = \\I - I, 



S\ 



1-1 



Proof: We have the following chain of equalities 

ml{\\A - P 5 /|| 2 : A G M g , trA = 1} = 
= M{\\A - Tf : A G M. q , tr A = 1 , T G T(S)} 

= inf | ^~J\F ■ A G M g , tr A ^ , T G T(g)| 

= m 4R™ :TGr( " 0) ' tr (/ TrfA )^ ^ 



sup 

5 50 



|tr(/TrfA)| 2 
f (TW 7 , J) 



{»^r (5 ») ltr (/T^o^ 



|P,o/||- 2 



□ 



For applications of the preceding theorems a description of the space A4(S°) d is needed. 
The identification of Ai(S°) d and of Ai(S c ) would be desirable from the point of view of 
prediction theory. Clearly, 



(3.9) M(S C ) C M(S°) d . 

However, whether equality really holds or not, seems to be a difficult problem, in general. 
It is related to basis properties of characters. Miamee and Pourahmadi ||MiP|| discussed 
this question and gave particular answers to it in case q = 1, cf. also JM3]. In contrast 
to these careful investigations, in the proofs of Thm. 4.1] and of ||CMiP| , Thm. 1] the 
equality Ai(S c ) = J\A(S°) d seems to be used in a rather general situation, but without 
any explanation how to prove it. 

We call a subset S C G to be G-exact if for any q G N and any W G >%(G*) the 
set identity JA(G\ S) = Ai(S) d holds. Because of the symmetry of this definition we 
immediately obtain the following fact. 

Lemma 3.8. The set S is G-exact if and only if the set G \ S is G-exact. 

The next result demonstrates that the alteration of a G-exact set S by finitely many 
elements does not have any influence on its G-exactness. 

Theorem 3.9. Let S be a subset of G and g G G \ S. Then S is G-exact if and only if 
S U {g} is G-exact. 
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Proof: Without loss of generality we may assume that G" S and g = 0. Suppose S 
to be G-exact. By flp) we have the chain of set inclusions M(S C ) C A^(5°) a! C A<(»S) a! C 
A4(G \ S 1 )- Thus, if G G .M(S ) d it can be decomposed as G = A + H for some A G M g 
and if G A4(S C ). So G - # = A G M(S°) d and hence, A = and G = H G A4(S C ). The 
latter means that S° is G-exact. To derive the G-exactness of S from the G-exactness of 
S° one has to combine the result just proved and Lemma [3.8| . □ 



Corollary 3.10. If either S or G \ S are finite subsets of G then S is G-exact. 

Proof: Define A4($) = {0}, so the empty set is G-exact by definition. With this 
setting the assertion follows from Theorem |3.9| and Lemma |3.8|. □ 



If G = Z then its dual group G* can be identified with the numerical interval [— 7T, n), 
where the addition of elements is understood to be done mod27r. The characters of 
G* = [— n, 7r) can be described as the set of functions e^(t) = e tkt , t G [— tt,it), k G Z. 
A univariate version of the assertion stated thereafter was given in [ MiP| , Lemma 3.6], 
compare with Thm. 3.1] for a univariate extension of this assertion to L p -spaces. 



Theorem 3.11. The set N is Z-exaci. 



Proof: By (|37 
M(Z\N). Then 

(3.10) 

and 

(3.11) 



we have A4(Z \ N) C A^(N) d . Let G G M(N) d be orthogonal to 



J e k G* dX = , fceN, 
J e k W~ 1 G*dX = , k G Z\N. 



By (|3.10| ) the element G* belongs to the Hardy space H 1 (of M g - valued functions), and 
( pj.ll[ ) implies that GW -1 belongs to H 1 , too. Therefore, GW _1 G* is an element of 
H 1 / 2 taking values in the set of non-negative Hermitian matrices. According to [ST] the 
function GW~ X G* is a constant function. However, by ( p,10| ) and ( |3.11| ) the index zero 
Fourier coefficient of GW^G* equals to zero, so GW^G* = and G = 0. □ 

4. Some applications 

From the general assertions of section three we can easily derive some prediction results 
for multivariate stationary sequences. Many of them are well-known, however here they 
are drawn from a different context. 



Theorem 4.1. (cf. JHe§ Thm. 4.3] or ]Ma|, Cor. 2.9]) 
Let W G >%(G*) and S = G \ {0}. Then 



PsI 

6 S 



W~ l d\ 
i 



W 



W~ x d\ 



and 



det ( / W- 1 dX 



-l/q 
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Proof: Since M(S°) d = {0} the first two results are immediate consequences of 



Theorem 3.2. The third result can be derived from 03 



□ 



Now let G = Z, and for a certain n G M set 

fi r 1 = Z\{0,l,2,...,n}. 



Theorem 4.2. (cf. |Ya|, Thm. 2] for the univariate case) 

Let n G N, W G W ? ([-7r,7r)) and Q = f e{W~ l d\, I G Z. Let the q x q-matrices 
{Dk : k — 1, n} be the solutions of the linear system 



(4.1) 
Then 



3 i 



n . 



k=l 



-1 r 



(4.2) 
and 



PsJ 



VSi 



V fc=i y 



J-^^efc)^" 1 



fe=i 



k=l 



'Si 



det / - ^ L> fc e fc , J 



fc=i 



-1/9 



Proof: According to Corollary |3.10| the set Si is Z-exact. Elementary calculations 
show that -Pgj/ = Sfc=i-Dfc e fc- Knowing this the result is an immediate consequence of 
Theorem EO and Lemma 13.41. □ 



Remark 4.3. From ( |4.2| ) we obtain 

n 

(4-3) E C- k Dl = C - A s l 



k=l 



Let T n be the Toeplitz matrix T n = (Cj-k)j,k=i,...,n, n G N. It is not hard to see that the 
system of equations ( |4.1| ) and ( |4.3| ) can be written in the form 

f —I \ ( -A Si x \ ( —I \ 



n+l ' 







i.e. 



( -A£ \ 




r- 1 . 

1 n+l 



VAW V o J \D* n j V o / 

Comparing the first q x g-blocks from both the sides we get / = OnA^ 1 , and therefore 
1 = det(fi n )det(A^ i 1 ). Here f2 n denotes the q x g-matrix in the left upper corner of the 
matrix r~ +1 . Using a well-known rule for computing minors of inverse matrices together 
with the fact that the nq x rag- matrix in the right lower corner of r n+ i equals T n , we get 

det(rg 



det(A 5l ) = det(O n ) 



det(r n+ i) 
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for the determinant of the prediction error matrix. 

For the special case where q = 1 Nakazi gave an elegant proof of Szego's infimum 
formula, cf. |N|, Cor. 4]. Zasukhin [g], Helson and Lowdenslager ||HelL| , Thm. 8] 
stated a version of Szego's theorem for non-negative Hermitian matrix-valued measures. 
Using the results of section three we can adapt Nakazi's proof to the multivariate situation 
to obtain the result by Zasukhin and Helson-Lowdenslager for absolutely continuous non- 
negative Hermitian matrix-valued measures. 

Theorem 4.4. Let G = Z and W be a function of L 1 , whose values are non-negative 
Hermitian matrices. Then 

{ 1 



(4.4) 



cxp 



log det WdX 



where the right-hand side of ( U-4\ ) has to be interpreted as zero in case the expression 
log det W is not integrable. 

Proof: First assume that W E W q ([-ir , it)) . If A E M' q and T E T(N) from 
the inequality between the arithmetic and the geometric means and from the fact that 
/log |det(A — T)\ d\ > log |det A\ = we obtain 



\A-T\ 



> 
> 



J tv[(A — T)W{A — T)*\ dX 
exp y logtrp - T)W(A - T)*] dX 
expjy log{detp-T)H/(A-r)*]} 1/(? dX 



Hence 
(4.5) 

Similarly, 
(4.6) 



> exp 



5 N > exp 



5^c > exp 



log det WdX 



log det W dX 



logdet^" 1 )^ 



By Theorem |3.11| and Theorem |3.9| the set N° is Z-exact, and by ( |3.7| ) the equality 
8j$ = S^c turns out. Comparing this with (|4.5|) and ( ^4.6|) we get the desired result under 
the additional assumption that W~ l E L . The general assertion can be derived from this 
partial one by a standard approximation argument demonstrated in [fA], pp. 260-261] for 
the case q = 1. □ 



From Theorem flj we obtain a somewhat surprising result. 
Corollary 4.5. Let G = Z and W E W 9 ([-tt, vr)). Then 

[det(A N )] lA? = inf y[det[(I -T)W(I -T)*]] 1/q dX : TeT(N)} 
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Proof: By and ((L4|) we conclude [det(A N )] 1/g = exp {l/q ■ J logdet W dX}. On 
the other hand, 

exp (l/q- J \ogdetWd\j = inf y [det[(I - T)W(I - T)*}} 1/q dX : T G T(N)| , 
cf. |K], Cor.]. □ 

5. Nakazi's prediction problem 

Another application of the results of Section three gives some information about the 
multivariate version of Nakazi's prediction problem. Let G = Z and S be of the form 

(5.1) S 2 = N c U{l,2,...,n} 

for some n G N. Assume the function W belongs to W g ([— 7T, 7r)) and is such that log det W 
is integrable. The set of such functions is denoted by VV'([— ir, n]) in the sequel. Let $ be 
the unique outer function of the Hardy space H 2 (of functions with values in M 9 ) such 
that 

(5.2) W = $*$ and J $ dX G M> . 

The function $~ x is also outer, but it does not belong to H 2 , in general. Let Bj be the 
j-th Taylor coefficient of j G N°. 



Lemma 5.1. (cf. ||MiP| , Cor. 3.7] for q 



Let W G W q ([—iT, 7r)). T/ien 

n n 

(5.3) P S2 = I- ^(BjB*)- 1 ■Y,n j < J i<\>- 

3=0 3=0 

and 



(5.4) 




Proof: By Lemma |3]^ and Theorems |3.9| and |3.11| the set S2 is Z-exact. Taking 
into account Theorem p.2| and the fact that the outer function belongs to H 2 as 
soon as W G W g ([— n, 7r)), the proof consists of straightforward calculations which will be 
omitted. □ 



Our goal is to establish (|5.3|) and (|5.4|) under the weaker assumption that logdet W is 
integrable. This can be done by an approximation procedure. Our approach is similar 
to that one presented in the proof of ||CMiP| , Th. 3] for the case q — 1. However, since 



we wish to compute not only the prediction error as done there but also the orthogonal 
projection, we give a complete proof of our generalized result. 

Let W m = W + 1/m ■ I, m G N. Denote by (., ,) m the M g -valued inner product 

of L 2 (W m ), by M m (S 2 ) the closure of T(S 2 ) in L 2 (W m ), and by P^' the orthogonal 
projection of L 2 (W m ) onto A4 m (S2). Note that L 2 (W m ) can be considered as a subset of 
L 2 (W), and that the inclusion 

(5.5) K(S 2 )CM(S 2 ), mGN, 
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Let $ m be the corresponding outer factor of W m , i.e 
W m = $l$ m and 
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lim (J - pjrh, I - P<g>I) m — (i — PsJ, i - PsJ) 



$ m d\ G M> , m G N . 



The following result is an easy consequence of [pLEfl , so the proof is omitted. 
Lemma 5.2. Let W G W q ([-n, it)). Then lim^oo $ m = $ w.r.t. the topology of L 2 . 
Let Bj m be the j-th Taylor coefficient of Set 



-i 



2_j B jm e j 
j=0 



and 



-i 



d=0 



3=0 



Lemma 5.3. Let W G W q ([—n,iv)). Then there exists a subsequence {nii}i eN ofN such 
that 



lim 

l— >oo 



p( m i) T 

^s 2 1 







Proof: Lemma [572] implies that there exists a subsequence {rrii}^^ of N such that 



lim 



$ a.e.. Since 



IWmX 1 ~ i^T'W 2 < J |$-$ m J 2 |$ m ;i 2 |$- 1 lW 1/ T^ = / \$-$ mi \ 2 \w-l\d\ 

and |$ - c&mj 2 !^ 7 ^ 1 ! < 2 ■ (|$| 2 + l^mj 2 )!^ 7 ""^ < 4 Lebesgue's dominated convergence 
theorem yields 



(5.7) 



lim II ($ 



1—fOO 



* \ — 1 
mi) 



0. 



Applying Lemma [54] we can write 

(5#-r($*)- 1 - = llT^y-i-T^r 1 ! 



T($ 



Lemma [5.2| implies that for every j G N° the j-th Taylor coefficient of $ m tends to the 
j-th Taylor coefficient of $ as m tends to infinity. Therefore, lim^oo T m{ = T uniformly 
on [— 7r,7r), and the second summand on the right hand side of ( |5.7| ) tends to zero as / 
tends to infinity. Also the first summand there tends to zero because of (|5.7| ) and the 
uniform boundedness of the set {T m; : I G N}. □ 



Theorem 5.4. Let G = Z, S 1 be the set described at (fjTjj and W G W'([— 7r, 7r)). T/ien 



-i 



(5.9) 



vi=o 
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(5.10) 

and 
(5.11) 



>s 2 



As 2 = ( E B ^ 



det (e^; 

u'=° 



-1/9 



Proof: From (|5.6|) and Lemma |5.3| we get 



I -P S2 I\\ > lim 11/ -P, 



( m 



2— »oo 



|/-(/-T($ 



If we take (|5.5| ) into account, we can conclude that (/ — T($*) r ) G A / f(S < 2). This yields 
( |5.9|) . To obtain ( |5.10| ) combine formula ( ^9|) with (^j]). Finally, (|5.11|) is a consequence 
of (p|). □ 

We conclude our paper by stating the prediction error matrix As for sets S of the form 

(5.12) 5 3 = N c U{n} 
and 

(5.13) S 4 = N c \{-n} 



for some n G N. The univariate versions of the results below can be found at [|CMiP| , 
Th. 5 and 6]. Let Aj be the j-th Taylor coefficient of the outer function $ which has the 
properties ( |5.2|) , j G N°. 

Theorem 5.5. Let G = Z, S 3 be the set described at ( \5.13j ) and W G W q ([-n, n)). Then 

-l 



(5.14) A 53 
// A n is regular then 



A*|/-A n [E^A] <\A . 



'n-l 



-1 



(5.15) 



A 53 = ^o(^)" 1 E A*Ai E A* A, A* n A Q . 



\j=o 



We omit the proof since ( |5.14|) can be obtained by a straightforward generalization of 
the proof given for ||CMiP| , Th. 5] to the multivariate case, and ( |5.15| ) follows by some 
simple matrix computations from ( |5.14|) . 



Theorem 5.6. Let G = Z, S A be the set described at ( \5.1Q ) and W G W'([-7r, 7r)). Then 



-1 



-1 



(5.16) 

// B n is regular then 
(5.17) 



A SA =A \l-B n \Y J B*B j ] B* n \ A*. 

Vi=o 



'n-l 



-1 



A^ 4 = MK)' 1 ( E B*B 3 E ) ■ 

3=0 I \i=o 
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Proof: Since the set 5*4 is Z-exact by Theorem |3 . 1 1| , Theorem |3.9| and Corollary p. 10 
and since the set S% is the reflection of 5*3 about the origin, the result immediately follows 



from the Theorems |3.2| and [5.5| whenever W E W 9 ([— 7r, tt)). If merely W E W'([— 7r, 7r)) 
approximate W by the sequence {W + 1/m ■ I} mt =n- 1=1 
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